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AIRCRAFT ENGINE FACE RADAR CROSS 
SECTION ANALYSIS 

BACKGROUND OF THE INVENTION 

The present invention is related generally to performing radar cross 
section (RCS) analysis of aircraft engine face components. Specifically, the present 
invention is related to a finite element based approach for the capturing of fine details 
of RCS patterns from large scale, complex engine face configxirations. 

5 The performance of a radar cross section (RCS) analysis on aircraft 

engine face components is a prohibitively difficult task due to the large amount of 
memory required to generate an appropriate electromagnetic description of the aircraft 
engine face components (problems) to be analyzed. The RCS analysis of aircraft 
engine face components usually requires special numerical and mathematical 

10 techniques to be developed to make the large problems practically solvable. These 
special numerical and mathematical techniques need to include appropriate handling 
of complex multi-layer design features, modeling of the entire axi-periodic structure 
of the aircraft engine face components and techniques for coupling a finite domain to 
an infinite domain. The coupling of the finite domain to the infinite domain is usually 

1 5 accomplished using an integral boimdary element technique. 

Numerous techniques for solving large-scale electromagnetic problems 
are in common use. These techniques include Method of Moments (an integral 
equation technique), finite element and finite difference techniques (differential 
techniques) and geometrical and physical theory of diffraction based techniques 

20 (asymptotic techniques). Attempts to use these traditional techniques to perform RCS 
analysis of aircraft engine face components have proven to be largely unsuccessful. 
The current techniques in use require large amounts of computer memory and 
computational time to generate results on simple models that do not correspond to 
actual devices or scenarios v^ith an aircraft engine face. The current techniques are 

25 essentially unworkable once they are applied to actual devices and scenarios because 
of the extraordinary amount of memory and computational time that is required to 
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model the spatial volume and resolve the complex material and geometric details of 
the device. 

Therefore what is needed in the art is a technique that can accurately 
perform an RCS analysis on actual aircraft engine face components v^thout requiring 
large and excessive amounts of computer memory and computational time. 

BRIEF SUMMARY OF THE INVENTION 

The present invention is directed to a computer program product, a 
system and a method of calculating, in the frequency domain, a radar cross section of 
an aircraft component having an axi-periodic structure. The present invention will 
first create a finite element model (commonly called a mesh or a grid) for a 
preselected period of an axi-periodic aircraft component. Next, a system matrix of the 
finite element model is transformed using a Discrete Fourier Transform into a 
plurality of independent modes in a block diagonal matrix. Each independent mode is 
then solved to determine that portion of the total electromagnetic near-field 
contributed by each mode. The solving of each of the independent modes involves the 
use of an "infinite" pipe to remove the effects of rim scattering. The electromagnetic 
far-field is calculated, using a well known integral equation, firom the electromagnetic 
near-field. The radar cross section for the aircraft component is then calculated, using 
a well known equation, from the electromagnetic far-field. Once the radar cross 
section for the aircraft component is determined, the aircraft component can be 
modified to achieve an optimal design. 

One advantage of the present invention is that it can accurately capture 
fine details of RCS pattems from large scale, complex engine face configurations. 
Another advantage of the present invention is that it is effective in the analysis of 
conductive, coated and composite material configurations and can be used for multi- 
configuration, multi-frequency and multi-polarization analysis of aircraft engine 
components. Still another advantage of the present invention is that the RCS analysis 
is solvable in practical central processing unit (CPU) time in a desktop workstation- 
type environment. Yet another advantage of the present invention is that the RCS 
analysis can be used by an aircraft engine face designer to achieve, in a reasonable 
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amount of time, a design with excellent RCS perfomiance. Since the techniques 
accurately correlate to measured RCS values, the amount of time and prototype costs 
required to optimize designs is substantially reduced. 

Other features and advantages of the present invention will be apparent 
from the following more detailed description of the preferred embodiment, taken in 
conjunction with the accompanying dravraigs which illustrate, by way of example, the 
principles of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 illustrates the basic design model used for the computer 
analysis of the RCS of aircraft engine face components. 

Figure 2 illustrates the axi-periodic structure of an aircraft engine face 

component. 

Figures 3A-3D illustrate schematically the steps involved in creating 
mathematical representations of pipes of preselected or "infinite" lengths. 

Figures 4A-4C illustrate schematically the steps used by the present 
invention to remove rim scattering effects. 

Figures 5 and 6 show a comparison between the RCS analysis of the 
present invention and the actual RCS measured for an aircraft engine face component. 

Whenever possible, the same reference numbers v^U be used 
throughout the figures to refer to the same parts. 

DETAILED DESCRIPTION OF THE INVENTION 

The present invention is directed to performing a computer simulation 
or analysis of the radar cross section (RCS) of aircraft engine face components. The 
computer simulation or analysis is performed in the frequency domain and can 
accurately model continuous wave behavior. The computer simulation or analysis can 
be executed on any general purpose computer having memory devices (RAM, ROM, 
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hard disk, CD-ROM, etc.), processing units (CPU, ALU, etc.) and input/output 
devices (monitor, keyboard, mouse, printer, etc.). In a preferred embodiment of the 
present invention, the computer simulation or analysis is executed on a "desktop" 
computer or workstation. Preferably, the "desktop" workstation has a multi-processor 
arrangement, wherein the number of processors preferably range from 5 to 20. The 
multi-processor arrangement can execute instructions either serially or in parallel, as 
required, to achieve the most efficient processing of instructions, hi another 
embodiment of the present invention, several "desktop" workstations or computers 
can be used in parallel to execute portions of the RCS computer simulation or 
analysis. 

Figure 1 illustrates the basic design model used for the computer 
simvilation or analysis of the RCS of preselected aircraft engine face components. The 
aircraft engine face components are those engine components that are first struck by 
radar and include vanes and their corresponding assemblies, fans and their 
corresponding assemblies, engine frames, compressors or any other aircraft engine 
component having an axi-periodic structure. While the present invention is directed 
to the computer analysis or simulation of aircraft engine face components, it is to be 
understood that the computer simulation or analysis of the present invention can be 
applied to any aircraft component that has an axi-periodic structure. 

The basic design model shovm in Figure 1 has a three-dimensional 
cavity 100, whose aperture 101 is in a particular plane in space 102 that intersects 
cavity 100. Within the cavity 100 is the aircraft engine face component or target (not 
shown), whose RCS is to be determined by the computer simulation. The computer 
simulation receives a mathematical representation of an incident wave 104 as an input 
at the aperture 101 of cavity 100. The incident wave 104 is a imiform plane wave or 
a transverse electromagnetic wave (TEM) impinged at a preselected angle and at a 
preselected frequency. Next, the computer simulation calculates the total 
electromagnetic field, the near-field, at the aperture 101 of the cavity 100. The total 
electromagnetic field is the sum of the incident wave 104 and the scattered waves 106 
from the rim represented by the intersection of plane 102 and aperture 101, the cavity 
100 and the preselected aircraft engine face component. The computer simulation 
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then determines the scattered electromagnetic field for the aircraft engine face 
component by removing from the total electromagnetic field, the incident wave 104 
and the scattered waves fi-om the rim and cavity 100. The scattered electromagnetic 
field fi-om the aircraft engine face component is then used as the basis for obtaining 
the RCS of the aircraft engine face component. 

To begin the computer simulation or analysis, a finite element mesh is 
created of the preselected aircraft engine face component to be analyzed. The aircraft 
engine face components to be analyzed will preferably have an axi-periodic structure. 
In other words, there is a period, slice or portion of the aircraft engine face component 
that is repeatable around a central axis to form the entire aircraft engine face 
component. To ftirther illustrate. Figure 2 shows an aircraft engine face component 
200 having eight vanes 202. The eight vanes 202 are spaced equally at 45 degree 
intervals about a central axis to cover an entire 360 degrees. This equal spacing of the 
eight vanes 202 permits any 45 degree slice 204 of the aircraft engine face component 
200 to be repeated or placed next to itself seven additional times to form the aircraft 
engine face component 200. It should be noted that the aircraft engine face 
component 200 illustrated in Figure 2 is meant to be illustrative and in no way 
limiting as to the type or structure of aircraft engine face component that can be 
analyzed by the present invention. 

Since the preselected aircraft engine face component to be analyzed is 
an axi-periodic structure, a discrete body of revolution theory can be applied to the 
aircraft engine face component so that only a period or slice of the aircraft engine face 
component has to be modeled with finite elements. By modeling only a period or 
slice of the aircraft engine face component, the amount of memory required for the 
finite element model is greatly reduced. In addition, calculations based on the finite 
element mesh or model of the slice are as accurate as calculations based on a finite 
element mesh or model of the entire aircraft engine face component. In other words, 
the application of the discrete body of revolution technique reduces the computational 
effort but does not decrease the accuracy of the model. 
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The finite element model or mesh of the slice of the preselected aircraft 
engine face component will use sufficient finite elements to accurately resolve the 
incident waves 104 that strike the model or mesh. The size of the wave from crest to 
crest is called the wavelength and is inversely proportional to the frequency of the 
wave. The number of elements used in creating the finite element model for a 
required frequency of incident wave 104 may often result in a finite element model 
with unknown variables numbering as large as one million or more. The size of a 
finite element used in creating the model is smaller than the wavelength of any 
incident wave 104 that is striking the model, and is sometimes much smaller than the 
wavelength in order to resolve fine geometric detail. In a preferred embodiment, the 
finite elements created and used are three to four times smaller than the wavelength of 
the incident wave 104 that is striking the model or mesh. The finite elements used to 
create the mesh or model of the slice of the preselected aircraft engine face component 
preferably have a brick, prism or tetrahedron shape, although other shapes can be used 
for the finite elements. In one embodiment of the present invention, the fmite 
elements used for the computer analysis are second order edge elements of the curl- 
conforming type. However, other orders (e.g. third order) and types of fmite elements 
can be used to create the model or mesh. Material properties of the aircraft engine 
face component such as its la (permeability) and e (permittivity) can be accounted for 
in the governing equations. In addition, the governing equations can also accoimt for 
the presence of R-cards or Z-cards that are thin sheets of material with conductive and 
dielectric properties that are placed on the aircraft engine face component. Finally, the 
use of finite elements permits the complex geometric structure of an aircraft engine 
face to be modeled with great fidelity. 

In the analysis, the Discrete Fourier Transform technique is applied to 
reduce the large, sparse system of equations resulting from the finite element assembly 
(discretization) to a block diagonal sparse system of equations. The matrix size for a 
single block of the block diagonal sparse system is 1/N of the original system, where 
N is the number of periods. The number of Discrete Fourier modes used in the 
transformation directly corresponds to the number of periods or slices required to 
generate the entire aircraft engine face component. For example, using the slice 204 
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shown in Figure 2, the slice 204 would have to be repeated or used an additional seven 
times for a total of eight to generate the entire preselected aircraft engine face 
component 200. Thus, applying a Discrete Fourier Transform to the system matrix for 
slice 204 would result in a block diagonal matrix having eight Discrete Fourier modes. 
Each Discrete Fourier mode in the block diagonal matrix represents a distinct portion 
of the electromagnetic field scattered by the entire preselected aircraft engine face 
component (360 degrees). Thus, the entire electromagnetic field in the cavity 100 is 
represented by the sum of the Discrete Fourier modes. Each Discrete Fourier mode is 
obtained from the solution of an independent block in the block diagonal matrix. 

As discussed above, each Discrete Fourier mode represents a portion of 
the total electromagnetic field throughout the cavity 100 from the preselected aircraft 
engine face component. However, for compactness and computational efficiency, a 
system matrix called an "impedance matrix" or "super-element" that represents the 
unknowns at the aperture 101 of the cavity 100 is all that is required for fiirther 
processing to obtain the RCS of the preselected aircraft engine face component. To 
obtain the impedance matrix at the aperture 101 of the cavity 100, each block is 
processed using a Back to Front solver. The Back to Front solver will process the 
unknowns in a block from the back (rear) of the cavity to the aperture 101 at the front 
of the cavity 100. During the progression from the aircraft engine face component to 
the aperture 101, the equations corresponding to unknovms in the intermediate layers 
are factored and then discarded because only the layer that corresponds to the aperture 
of the cavity 100 is necessary. At the end of this process, ah impedance matrix for the 
unknowns on the aperture 101 of the cavity 100 is obtained. This impedance matrix 
accurately represents the response of the entire cavity 100 and the preselected aircraft 
engine face component to an impinging electromagnetic field or incident wave. One 
impedance matrix is obtained for each Discrete Fourier mode. The impedance matrix 
can also be called a "super-element" because it represents a finite element model of 
the entire cavity 100 and the preselected aircraft engine face component which can be 
repeatedly used for more than one analysis. This super-element is stored in computer 
memory or on hard disk drive storage for fiirther processing to be described shortly. 
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To accurately represent the electromagnetic fields observed in the test 
environment, two additional tasks must be completed: the modeling of a zero duct 
length (ZDL) test fixture and the removal of the effects of the rim. These two 
additional tasks will be later combined v^th the cavity impedance matrix generated 
above to analyze the entire electromagnetic system. 

The ZDL test fixture is an actual physical device in which the 
preselected aircraft engine face component is placed for measurement of the RCS 
jfrom the aircraft engine face component. An impedance matrix for the ZDL test 
fixture has to be created for each Discrete Fourier mode of the computer simulation or 
analysis so that the RCS determined by the computer simulation corresponds to the 
actual test setup. 

The steps in creating an impedance matrix that represents the ZDL test 
fixture are illustrated in Figures 3A-3D. As shown in Figure 3A, the first step is to 
create a fmite element model 300 of a preselected sublength L of the ZDL fixture and 
generate its resulting system matrix. Next, as illustrated in Figure 3B, the ZDL system 
matrix is reduced into an impedance matrix 302 by eliminating the unknowns in the 
interior. The resulting impedance matrix 302 contains the unknovms on the two end 
surfaces of a ZDL test fixture of preselected sublength L. This impedance matrix 302 
also includes the coupling information between unknowns on the two end surfaces. 
The process is continued recursively as shovm in Figures 3C and 3D to obtain the 
impedance matrix for any fmite length of the ZDL test fixture. In Figure 3C, an 
impedance matrix 304 for the ZDL test fixture of length 2L is obtained by recursively 
adding two impedance matrices for the ZDL test fixture of length L and then reducing 
the resulting matrix to eliminate the interior unknovms. The process is then repeated, 
as shown in Figure 3D, where the impedance matrix 304 is added to itself and then 
reduced into an impedance matrix 306 representing the two end surfaces of a ZDL test 
fixture of length 4L and the coupling information between unknowns on the two end 
surfaces. This process is repeated until the preselected length of the ZDL test fixture 
is achieved. Once the preselected length of the ZDL test fixture is achieved, this ZDL 
impedance matrix, or "super-element" is stored in computer memory or on the hard 
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drive for later use. One impedance matrix for the ZDL test fixture is obtained for each 
Discrete Fourier mode. 

An important aspect of the present invention is the construction and 
use of a "infinite pipe" to represent the electromagnetic fields extending from the 
cavity 100 containing the preselected aircraft engine face component and ZDL test 
fixture to very far distances. Without this infinite pipe representation, the fields 
would "bounce" or reflect back from artificially imposed structures in the numerical 
model, thus corrupting the solution and reducing the accuracy of the RCS 
computation. In other words, the infinite pipe eliminates the reflections from 
artificially imposed structures. 

The infinite pipe is a numerical model of a very long pipe or cavity 
filled with an extremely low loss material to dampen the propagating electromagnetic 
field, thereby eliminating any reflection from the far end of the pipe. The infinite pipe 
starts with a system matrix from a finite element model of finite length, as was done 
with the ZDL model. Often, the ZDL super-element is used to model this first stage 
of the infinite pipe, although this is not necessary. The length of the infinite pipe is 
then recursively enlarged using the process described with regard to Figures 3C and 
3D. As opposed to the process for creating the ZDL test fixture impedance matrix, the 
infinite pipe process is continued until the coupling terms between the unknowns on 
opposite surfaces of the infinite pipe impedance matrix or super-element are 
negligible. This process causes the opposite surfaces in the infinite pipe impedance 
matrix to become uncoupled. Only a few iterations are required because the length of 
the pipe grows rapidly as a function of 2", where n is the recursive iteration number. 
Once completed, the terms containing unknowns on the far-side of the infinite pipe 
super-element are discarded. The near-surface terms of the infinite pipe super- 
element, i.e. the impedance matrix adjacent to the front surface of the ZDL, are then 
stored in computer memory or hard disk storage for later combination with the ZDL 
and cavity super-elements. Again, one infmite pipe impedance matrix or super- 
element is stored for each Discrete Fourier mode. 



13DV-13689 




The impedance matrices or super-elements for the ZDL test fixture, the 
infinite pipe and the cavity are further processed to obtain the RCS of the aircraft 
engine face component. The processing of these super-elements or impedance 
matrices involves the coupling, complete elimination, factorization and solving of the 
matrices in a process to be described below. The construction of separate super- 
elements for different portions of the computer analysis or simulation process permits 
possible "swapping" of super-elements for an iterative design analysis. That is, 
storing the super-elements for future use can save computational time and effort when 
applied to other configurations. 

The process for calculating the RCS of an aircraft engine face 
component begins by determining the total electromagnetic field at the interface of the 
ZDL test fixture and infinite pipe. Figures 4A-4C illustrate schematically the process 
for determining the scattered electromagnetic field for an aircraft engine face 
component and for removing the effects of "rim" scattering. Figure 4A illustrates the 
configuration for obtaining the total electromagnetic field for each Discrete Fourier 
mode. In Figure 4A, the infinite pipe 400 is combined at interface or junction 402 
with the cavity/aircraft engine face component 404 and the ZDL test fixture 406. The 
impedance matrices for the cavity/aircrafl engine face component 404 and the ZDL 
test fixture 406 are coupled and the interior unknowns eliminated to generate a new 
impedance matrix representing both the cavity/aircraft engine face component 404 and 
the ZDL test fixture 406. The new unpedance matrix representing both the 
cavity/aircraft engine face component 404 and the ZDL test fixture 406 is then 
coupled to the impedance matrix of the infinite pipe 400 to generate a final impedance 
matrix. The fmal impedance or system matrix for the configuration of Figure 4A is 
then factored and solved using an incident plane wave of desired fi-equency and 
direction impinging at the interface 402 between the ZDL test fixture 406 and infinite 
pipe 400 and returns the total electromagnetic field (near-field) contributed by that 
Discrete Fourier mode. The scattered electromagnetic field (near-field) contributed by 
a Discrete Fourier mode for the akcraft engine face component is then determined by 
subtracting the incident field from the calculated total field at interface 402. The 
scattered electromagnetic far-field contributed by a Discrete Fourier mode for the 
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aircraft engine face component is then determined using a well known integral 
equation from the calculated near-field. The electromagnetic far-field is the 
electromagnetic field present at a substantial distance from the electromagnetic near- 
field. Finally, the RCS is determined using a well known equation from the sum of 
the scattered electromagnetic far-fields contributed by each Discrete Fourier mode. 
However, in using the configuration shown in Figure 4A there is an electromagnetic 
field generated from the scattering of the incident wave from the rim or boundary 
between the ZDL test fixture 406 and the infinite pipe 400. To obtain an accurate 
computer analysis of the RCS of an aircraft engine face component, this rim scattering 
effect must be removed. The infinite pipe super element is used to remove the effects 
of rim scattering in a process described below. 

As shown in Figure 4B, a similar configuration is generated to the one 
in Figure 4A. However, in Figure 4B, the infinite pipe super-element is coupled to an 
identical infinite pipe super-element and processed as described above to obtain a 
final impedance matrix for the configuration. The second infinite pipe 408 represents 
a domain extending a large distance opposite the original infinite pipe 400 and 
replaces the ZDL test fixture 406 and cavity/aircraft engine face component 404 
sections shown in Figure 4A. The same incident plane wave used to solve the fmal 
system matrix for the configuration in Figure 4A is impinged at the intersection 410 of 
the two infinite pipes. The final impedance matrix for the system of two-infinite pipes 
is solved, as described above, for the electromagnetic fields (near-fields) at 410 and 
then the electromagnetic far-fields are determined as described above. The rim 
scattering effect is then reduced to sufficiently low values by subtracting the far-fields 
from the two infinite pipe solution of Figure 4B from the far-fields solution of Figure 
4A. This difference yields the scattered electromagnetic field for the aircraft engine 
face component 404 as shown in Figure 4C. The RCS is then calculated from the 
remainder of these two subtracted solutions. 

However, in another embodiment of the present invention, the rim 
scattering effect can be removed and the scattered electromagnetic field for the aircraft 
engine face component can be obtained by subtracting the electromagnetic near-field 
solution for the two infinite pipe solution of Figure 4B from the electromagnetic near- 
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field solution for the configuration of Figure 4A. The scattered electromagnetic field 
(near-field) for each of the independent Discrete Fourier modes can then be added to 
determine the total scattered electromagentic near field. The total scattered 
electromagentic far-field is determined firom the total scattered electromagnetic near- 
5 field. Finally, the RCS is determined firom the total scattered electromagnetic far- 
field. 

In another embodiment of the present invention, the use of conjugate 
mode symmetry can reduce the computational processing of the Discrete Fourier 
modes. Conjugate mode symmetry involves the use of primary modes and conjugate 
modes that are directly related to the primary modes. Conjugate mode symmetry 
permits the impedance matrices for certain modes to be reused for other modes. The 
impedance matrix for a conjugate mode is the transpose of the impedance matrix for 
the primary mode. Conjugate mode symmetry computations can reduce the number of 
required mode impedance matrices to (N+2)/2 (for odd values of N, eliminate the Vi 
remainder), where N is the original number of modes. For example, assume that an 
aircraft engine face component would generate five Discrete Fourier modes. Applying 
the conjugate mode theory to these five modes would result in the first mode being 
distinct, the fifth mode being a conjugate of the second mode and the fourth mode 
being a conjugate of the third mode, thereby reducing the total number of modes to be 
independently computationally processed to three. Thus, by using the conjugate mode 
theory not all of the modes have to be independently processed thereby significantly 
reducing the amount of computational time required to generate a solution. 

After the RCS is calculated for the aircraft engine face component, 
changes can be made to the design of the aircraft engine face component to obtain an 
25 optimal RCS for the aircraft engine face component. Some of these changes can 
include a slight alteration of the material properties (e.g. ^i, e), the addition or removal 
of material treatments or a slight modification to the geometry of the aircraft engine 
face component. The computer simulation can then be executed' again to obtain a new 
RCS value for the modified aircraft engine face component. This modification and re- 
30 execution process can be repeated until the optimal RCS value for an aircraft engine 
face component is obtained. Each new design configuration of the aircraft engine face 
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component does not require new infinite pipe and ZDL test fixture super-elements to 
be created unless the number of Discrete Fourier modes is changed. 

Figures 5 and 6 show a comparison of the RCS determined by the 
computer simulation of the present invention and the RCS measured using a ZDL test 
fixture for several different frequencies and polarizations. Figure 5 show a 
comparison of measured and analyzed RCS values for different frequencies of 
incident waves. The incident waves used for the comparisons in Figure 5 have an HH 
(horizontal) polarization and the ZDL test fixture has a thirty inch diameter. Figure 6 
also shows a comparison of measured and analyzed RCS values for different 
frequencies of incident waves. The incident waves used for the comparisons in Figure 
6 have a W (vertical) polarization and the ZDL test fixture has a thirty inch diameter. 
As can be seen from Figures 5 and 6, the computer simulation of the present invention 
can accurately simulate what RCS will be measured for a particular aircraft engine 
face component in the ZDL test fixture. 

While the invention has been described wdth reference to a preferred 
embodiment, it will be understood by those skilled in the art that various changes may 
be made and equivalents may be substituted for elements thereof without departing 
from the scope of the invention. In addition, many modifications may be made to 
adapt a particular situation or material to the teachings of the invention without 
departing from the essential scope thereof Therefore, it is intended that the invention 
not be limited to the particular embodiment disclosed as the best mode contemplated 
for carrying out this invention, but that the invention will include all embodiments 
falling within the scope of the appended claims. 
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